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The  response  of  titanium  alloy  (6A1-4V)  to  simulated  nuclear 
cloud  particle  environments  is  addressed  in  this  combined  experi- 
mental and  analytical  study.  The  results  are  applicable  for 
ascent  vehicles  with  hot  titanium  structures  flying  through 
nuclear  burst  ice  or  dust  clouds.  A three-phase  exploratory 
program  was  conducted:  An  initial  oxidation  test  series  using 
a»  SAI-designed  test  fixture,  a test  series  using  the  AEDC  Dust r j 


DD  1473^  EDITION  OF  I NOV  «5  IS  OBSOLETE 


security  classification  of  this  page  Dara  Enfaratf) 


H'OI 


aA'  r 


SECURITY  CLASSIFICATION  OF  THIS  PAGEflWijn  Dmim  Enlarmd) 
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Erosion  Tunnel  (DET) , and  a particle  impact  test  series  using 
the  SAI  hypervelocity  impact  facility.  Analytical  modeling 
of  titanium  oxidation  produced  results  in  close  agreement  with 
experiment.  Although  increased  oxidation  of  titanium  was  con- 
firmed in  HVO  vapor  versus  oi'  environment  tests,  no  net  effect 
of  liquid  H5O  and  0^  oxidation  under  combined  particle  erosion 
conditions  was  observed  in  the  AEDC  DET.  This  is  believed  to  be 
the  result  of  the^heat  of  vaporization  of  H^O  exceeding  the  heat 
of  oxidation  of  H^O  vapor. 
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INTRODUCTION  AND  SUMMARY 


Extensive  testing  and  analyses  (1-5)*  have  been 
conducted  to  develop  a response  model  (erosion  and  heating) 
for  titanium  subjected  to  dust  impacts.  This  model  has  been 
successful  in  a broad  range  of  ground  test  facilities  in 
predicting  titanium  response.  However,  consideration  has 
not  been  given  to  the  combined  effects  of  nuclear  cloud 
ice/water  particle  and  dust  particle  impacts  on  titanium. 

In  general,  dust  particle  impacts  would  be  expected  to  be  an 
over-simulation  of  ice  impacts  and  thus,  a design  to  a dust 
environment  of  equivalent  mass  concentration  would  be  con- 
servative. However,  there  remains  some  uncertainty  relative 
to  the  effect  of  ice/water  impacts  on  titanium  surface 
chemistry  heating  augmentation.  This  document  describes  a 
combined  analytical  and  experimental  program  to  bound  the 
effects  of  ice/water  on  titanium  surface  chemical  reaction 
heating  augmentation  and  to  develop  an  improved  heat  transfer 
model  to  describe  these  effects,  if  necessary,  for  design 
applications. 

1.1  Technical  Background 

Extensive  studies  of  titanium  oxidation  have  shown 
the  dominant  surface  chemistry  heating  augmentation  mech- 
anism to  be  the  result  of  the  formation  of  solid  state 
oxides  (3,  4).  Oxygen  available  at  the  surface  of  the 
titanium  diffuses  through  a protective  oxide  coating  and 
releases  energy  in  the  form  of  a heat  of  solution.  As  the 
temperature  of  the  titanium  or  the  oxygen  concentration 
increases,  the  oxygen  diffuses  faster  and  energy  is  released 
at  an  increased  rate.  Titanium  ignition  occurs  when  the 
heat  released  due  to  the  chemical  reaction  exceeds  the  heat 
losses  due  to  convection,  conduction,  and  radiation.  Titanium 
nosetips  must  be  designed  to  stay  well  below  this  ignition 


* Refers  to  reference  numbers  in  Section  6.0. 
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* temperature  threshold.  Particle  impacts  complicate  the 
oxidation/heat  release  process  since  they  remove  the  pro- 
tective oxide  coating  and  increase  the  temperature  locally, 
thus  increasing  the  local  rate  of  oxygen  diffusion  and 
chemical  heat  release.  Figure  1-1  describes  this  process 
schematically. 

Past  studies  (4)  have  produced  models  for  the  chem- 
ical heat  release  in  titanium.  Figure  1-2  shows  schemati- 
cally that  the  rate  of  oxygen  diffusion  into  titanium  (and 
thus  rate  of  chemical  heat  release  q . ) is  controlled  by 
the  mass  flux  of  oxygen  at  the  surface  (m  ) (in  the 

boundary  layer  flow  field)  and  by  the  rate'^of  oxygen  dif- 
fusion in  the  solid  titanium  (m^  ) dominating 

rate  process  is  a function  of  the  time  between  impacts  t 
(obscuration  time) . The  chemical  reaction  heat  release  can 
be  written: 


where  AH^  is  the  heat  of  solution.  For  typical  dust  cloud 
conditions,  the  process  is  dominated  by  solid  phase  dif- 
fusion (this  must  be  reexamined  for  nuclear  ice  cloud  con- 
ditions) . The  oxygen  diffusion  into  the  solid  can  be  written: 


where  A and  E are  constants  related  to  diffusion  of  oxygen 
in  titanium,  T is  the  titanium  temperature,  t is  the  time, 
and  Cg  is  the  oxygen  concentration  at  the  surface.  Com- 
bining the  above  expressions  yields  the  functional  form 
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FIGURE  1-1.  SCHEMATIC  OF  TITANIUM  SURFACE  REACTION  PROCESS 


FIGURE  1-2.  SCHEMATIC  OF  CONSTRAINTS  ON  OXYGEN  MASS  FLUX  INTO 
SOLID  TITANIUM 


/? 


exp  (-B/T) 


’react  ' 


2 p d 


where  B is  a constant  and  t = 


where  a is  the  ratio 


^ xo  o uciii  u cliita  «.  — ^ wiicxc  ix  xo  UAic  xauxw 

3 Priy  « 

of  projected  crater  surface  area  to  particle  cross-sectional 
area  (order  of  1),  is  the  particle  density,  d^  is  the 
diameter,  V is  velocity  and  p^j^  is  the  cloud  concentration. 
This  model  is  generally  applicable  for  impacts  of  any  part- 
icle type  for  oxygen  diffusion  into  titanium  from  any  oxygen 
source.  The  constants  in  the  above  relation  must  be  deter- 
mined for  specific  cases.  For  dust  impacts  on  titanium  in 
air  environments,  a suitable  empirical  relation  was  derived 
to  be 


5 , . = 47700  — exp  (-23085/T (°R) ) Btu/ft^-sec 

r GaC  u 

where  T is  assumed  to  be  the  local  titanium  surface  temper- 
ature. A comparison  of  past  AEDC  DET  data  (4)  with  this 
model,  shown  in  Figure  1-3,  shows  excellent  agreement.  The 
more  generalized  model  for  dust  particle-titanium  chemistry 
coupling  provides  a framework  for  considering  the  additional 
effects  of  ice  impacts  on  titanium  surface  chemistry. 

There  are  several  unique  differences  between  ice  and 
dust  impact  conditions  which  could  alter  the  titanium  sur- 
face chemistry  heating  augmentation.  Specific  parameters 
affected  include  t (particle  density,  diameter  and  cloud 
density  differences) ; T (different  kinetic  energy  transfer 
and  thus  local  crater  temperature  on  impact) ; (available 
concentration  of  oxygen  in  the  boundary  layer  at  the  titanium 
surface) ; B (activation  energy)  and  the  rate  constant  for 
diffusion  of  water  species  into  titanium. 


The  objective  in  this  program  was  to  obtain 


oxidation  and  impact  test  data  for  titanium  at  elevated 
temperatures  in  oxidizing  environments  to  determine  if  the 
presence  of  ice  or  water  significantly  altered  the  response 
characteristics  of  titanium  under  dust  impact  conditions 
which  were  successfully  modeled  in  previous  studies.  Test 
facilities  selected  were  as  follows: 

1.  SAI  Oxidation  Test  Facility 

2.  SAI  Ice  Impact  Test  Facility 

3.  AEDC  Dust  Erosion  Tunnel 

No  facility  currently  exists  v/hich  can  simulate  multiple  ice 
erosion  on  titanium  in  an  oxidizing  environment  for  time  periods 
greater  than  a few  milliseconds  in  order  to  evaluate  chemistry 
effects.  The  above  facilities  simulate  critical  elements  of  the 
nuclear  cloud  ice/water  and  dust  environment  for  evaluating 
a missile  titanium  component  response.  In  addition  to  the 
test  program,  an  improved  analytical  model  was  formulated 
for  titanium  oxidation  based  on  data  in  the  literature 
including  limited  data  on  the  effect  of  steam  oxidation  of 
titanium.  This  model  was  used  to  analyze  test  data  and 
evaluate  the  combined  effects  of  ice/water  and  dust  particle 
environments  on  titanium. 

1.2  Summary  of  Results 

The  results  of  this  exploratory  program  show  no 
significant  overall  degradation  in  titanium  performance  for 
combined  ice/water  and  dust  environments  versus  dust  en- 
vironments. Previous  studies  of  titanium  response  to  dust 
showed  surface  chemical  reactions  becoming  important  at 
temperatures  in  excess  of  2000°F  and  eventually  causing 
ignition  at  higher  temperatures.  This  temperature  of  2000®F 
remains  a valid  criteria  for  ice/water  and  dust  below  which 
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surface  chemistry  effects  are  unimportant.  Ignition 
is  not  expected  where  the  bulk  temperature  of  the 
titanium  is  always  below  2000®F. 

Specific  results  from  this  program  contributing  to 
the  above  general  conclusion  include: 

1 . Titanium  Oxidation  Results 

The  SAI  Phase  I tests  of  titanium  in  non-erosive 
oxidizing  environments  (H2O  vapor  and  ©2)  showed  an  apparent 
increase  in  oxidation  rate  in  H2O  vapor  compared  to  for 
temperatures  of  1200-1500®K.‘  At  the  1700°K  condition,  this 
trend  was  not  evident  with  O2  and  H2O  vapor  giving  similar 
oxidation  rates.  Figure  1-4  shows  a comparison  of  the  Phase 
I experimental  results  with  the  oxidation  model  derived  from 
the  literature.  The  nominal  model  suggests  a parabolic-type 
weight  gain  law  (AW//t)  where  AVJ  is  the  total  weight  gain 
and  t is  the  total  test  time.  The  bulk  of  the  test  data 
(both  H2O  vapor  and  ©2)  fall  within  the  + factor  of  2 un- 
certainty band  for  the  ©2  oxidation  model  based  on  data  in 
the  literature.  The  exception  is  the  1.0  atm  data  at  1400°K 
where  a pressure  effect  is  suggested.  This  was  unexpected 
based  on  a solid-phase  diffusion-controlled  model  and  sug- 
gests the  possibility  that  the  SAI  oxidation  tests  may  have 
been  gas-phase  rate-limited  although  this  was  not  confirmed 
during  the  testing.  Although  the  H2O  vapor  data  from  this 
program  falls  within  the  uncertainty  band  of  the  total 
experimental  data  base  for  O2  oxidation  data  from  the  litera- 
ture, the  current  K2O  vapor  data  does  indicate  a general 
increase  in  oxidation  compared  to  the  current  ©2  data  up  to 
the  1700°K  condition. 
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FIGURE  1-4.  COMPARISON  OF  EXPERIMENTAL  AND  ANALYTICAL 
TITANIUM  OXIDATION  RESULTS 
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2.  Combined  Erosion  and  Oxidation 


The  Phase  II  AEDC  Dust  Erosion  Tunnel  (DET)  test 
results  show  no  net  increase  in  heat  transfer  for  titanium 
models  in  combined  water  vapor,  air  and  dust  environments 
versus  air  and  dust  environments.  Figure  1-5  illustrates 
this  fact  for  a typical  DET  run.  The  backface  titanium 
hemisphere  temperature  versus  time  is  shown  for  four  models 
in  different  environments  --  clear  air,  clear  air  plus  H2O, 
dust,  and  dust  plus  H2O.  The  effect  of  the  H2O  environment 
in  all  cases  was  a net  cooling  effect.  Although  some  in- 
creased oxidation  effect  due  to  H2O  vapor  may  be  present  (as 
indicated  above) , the  heat  of  vaporization  dominates  this 
oxidation  heat  release  resulting  in  the  overall  cooling 
effect.  This  has  been  confirmed  in  other  DET  runs  up  to 
bulk  titanium  temperatures  of  2500“F. 

3 . Ice  Impacts  in  Oxidizing  Environments 

The  Phase  III  SAI  results  compared  single  ice  and 
glass  impacts  on  titanium  under  oxidizing  conditions  at  a 
bulk  temperature  of  'V1700°F  ('^1200°F)  . Figure  1-6  illus- 
trates the  surface  characteristics  of  two  typical  test 
specimens.  Although  the  ice  produced  less  damage  (smaller 
crater)  than  the  glass  particle,  both  particle  types  pro- 
duced spallation  of  the  protective  oxide  coating.  However, 
no  ignition  was  observed. 


1-10 


ru:j 


w 

ft,  c:  o 
D M tN  o 
o <c  rn  <N 

u 'V  i:: 

O CC  Du  \ 

iu  < < H t" 

cr!  U tJ  to  00 

U J o a 

=:  u u a Q 


0 

4J 

U 

0 

CO 

0 

m 

\ 

to 

II 

m 

G 

0 

u* 

ca 

cu 

IT) 

Q 

II 

0 

u 

0 

0 

(U 

0 

CO 

ri 

n! 

tn 

m 

1 

•H 

\ 

\ 

0 

u 

c: 

0) 

XI 

XI 

rH 

u 

Q- 

Cu 

rH 

c« 

0 

N 

u 

m 

cx> 

in 

CD 

0 

V 

<N 

VO 

0 

m 

• 

< 

• 

in 

• 

0 

<T\ 

M 

m 

rH 

rg 

• 

m 

Q 

Q 

U 

II 

II 

II 

II 

II 

h-i 

rn 

0 

U 

0 

[-• 

M 

< 

< 

o: 

0 

CO 

Q 

M 

[-• 

£-• 

M 

(N  0 

0 

1 

tn 

cn 

< 

Q 

<N 

cu 

C- 

Eh  : : = : 

rr  O <Ti  CT\ 
M ^ in  *3*  rr 

o o O O O 


(N  ID 

.-I  r-'  <-l  rH 

I I I I 

-»H  H -H  "iH 

Eh  E-<  Eh  H 


&;OD«> 


< s 

o 7:: 
a o 

< E- 

u 00  o 
u Q r: 


••  D •• 

” M o 
U ^4  < 
CU  O iJ 
O to  Cm 


in 

U 

2 

u 

D 

to 

q:: 

Eh 

*-F 
« -< 

U 

VO 

0 

M 

Q 

Eh 

U 

u 

0 

w 

U 

< 

M 

s 

in 

0 

cx: 

tJ 

Q 

Ua 

g 

CO 

Eh 

7^ 

D 

d 

CO 

Ci 

u 

lu 

a: 

H 

h-t 

in 

ro 

r'j 

tH 

FIGURE  1- 

Jo  Juti.LVUJJKJ.i.  JOVJJ.'ivsi 


2.0 


SAI  OXIDATION  AND  IMPACT  TEST  SERIES 


2.1  Test  Objectives  and  Scope 

The  objectives  of  the  Phase  I oxidation  tests  were  to 
evaluate  the  effect  of  water  vapor  and  pressure  on  the 
oxidation  characteristics  of  titanium  (6AL-4V)  at  elevated 
temperatures.  The  scope  of  Phase  I included  three  test 
series  with  a total  of  36  tests  using  the  oxidation  test 
fixture  at  the  SAI  Santa  Ana,  California  laboratory. 

The  objective  of  the  Phase  III  impact  tests  was  to 
evaluate  titanium  response  to  ice  and  glass  impacts  in 
oxidizing  environments  to  determine  if  ignition  occurs 
from  single  impacts.  The  scope  of  Phase  III  included 
six  tests  in  the  SAI  impact  facility  at  Santa  Ana,  Califor- 
nia on  6AL-4V  titanium  at  elevated  temperature  comparing  glass 
and  ice  impact  results. 

2.2  Oxidation  Test  Program 

The  Phase  I oxidation  tests  were  conducted  at  the  SAI 
Santa  Ana,  California  laboratory  in  a test  rig  especially 
designed  for  this  DNA  program.  The  test  rig  is  shown  in 
Figure  2-1.  The  test  specimen  (see  Figure  2-2)  is  suspended 
by  thermocouple  wires  welded  to  the  specimen  in  a vicor 
cylinder.  An  induction  coil  attached  to  the  25KW  induction 
heater  surrounds  the  vicor  tube.  The  environment  (O2/  air, 
argon,  H2O  vapoi , or  mixtures)  is  controlled  by  admitting 
the  respective  gases  to  the  test  chamber  via  a control  panel 
equipped  with  solenoid  on-off  valves  and  manually  operated 
flow  regulating  valves. 

The  insitu  instrumentation  consisted  of  type  E and 
type  S thermocouples  and  an  optical  pyrometer  to  measure 
specimen  temperature  history.  Pre-test  and  post-test  weight 
measurements  of  the  test  sample  and  its  thermocouple  were 
performed  to  determine  oxidation  weight  gain. 
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FIGURE  2-lb.  OXIDATION  TEST  STAND  - REAR  VIEW 


The  test  specimen  configurations  are  shown  in  Figure 
2-2.  The  disc  configuration  with  a Type  E thermocouple 
welded  to  the  backface  was  used  for  check-out  and  the  Phase 
III  impact  t s.  The  rectangular  specimens  with  a Type  S 
thermocouple  welded  to  the  end  were  used  exclusively  for  the 
oxygen  weight  gain  tests. 

The  test  procedure  followed  to  obtain  approximately 
isothermal  weight  gain  data  consisted  of  the  following 
steps;  1)  a pre-test  weighing  of  the  test  specimen;  2) 
installing  the  test  specimen  and  instrumentation;  3)  es- 
tablishing the  desired  pressure  of  the  oxidizing  environment 
with  a continuous  flow  of  the  same  gas  through  a test  fix- 
ture; 4)  initiating  the  induction  heater  and  adjusting  the 
power  to  attain  the  desired  temperature  level  based  on  real- 
time observation  of  the  specimen  temperature  history;  5) 
maintaining  constant  temperature  for  the  prescribed  test 
time;  6)  cut-off  of  oxidizing  environment  and  power;  and,  7) 
post-test  weighing  of  the  test  specimen  plus  any  oxide  spall 
(captured  in  a catcher  at  the  bottom  of  the  test  rig) . 

Figure  2-3  shows  thermocouple  and  pyrometer  test 
results  for  a typical  low  and  high  temperature  test  run. 

The  thermocouple  measurements  constitute  the  primary  data 
with  the  pyrometer  as  back-up.  Although  a comparison  of  the 
temperatures  indicated  by  the  thermocouple  and  the  pyrometer 
is  not  possible  without  emissivity  data,  confidence  in  the 
data  can  be  gained  by  comparing  the  temperature  changes 
indicated  by  the  two  independent  instruments  during  the 
"constant"  temperature  period  of  the  test.  To  this  end,  the 
relative  sensitivities  of  the  two  instruments  are  indicated. 

The  results  of  the  Phase  I oxidation  tests  are  tabu- 
lated in  Table  2-1  for  three  series  of  tests.  Appendix  B 
provides  detailed  tabulated  temperature-time  data. 
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NOTE:  All  specimens  1 cm  x 2 cm  x 0 . 2 cm,  Mass  = 1.8  gms. 

Reference  Temperature,  18° C.  Chart  speed  0.1  in/sec 
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The  interpretation  of  the  data  and  the  identification 
of  specific  trends  is  deferred  to  Section  4 where  the  data 
is  presented  in  light  of  predictive  oxidation  models.  Two 
general  observations  which  can  be  made  from  the  raw  data  are 
depicted  in  Figure  2-4  where  results  obtained  at  1500°K  are 
presented.  A significant  increase  in  oxygen  uptake  is  seen 
for  tests  in  water  over  tests  in  pure  oxygen  at  1500°Kelvin. 

A similar  trend  was  observed  at  1200°K,  however,  no  difference 
was  measured  at  1700°K.  It  is  also  surprising  to  note  that 
data  obtained  with  a 50  percent  mixture  of  oxygen  and  water 
vapor  agreed  with  data  obtained  in  pure  oxygen  at  all  tempera- 
tures tested.  Unfortunately,  the  limited  data  generated  does 
not  permit  any  conclusions  to  be  drawn  from  these  observations. 

For  the  two  pressure  levels  tested,  namely  1/2  and  1 
atmosphere,  both  water  and  oxygen  data  suggest  an  increase  in 
oxygen  uptake  with  increasing  pressure.  This  trend  was 
unexpected  although  reported  once  in  Ref.  10  since  gas  phase 
diffusion  cannot  be  important  in  a single  component  atmosphere. 
As  will  be  discussed  in  Section  4.0,  this  pressure  effect 
may  result  from  gas  phase  rate-limiting  processes  as  opposed 
to  any  basic  change  in  the  solid  phase  oxidation  or  dissolu- 
tion processes. 

During  the  Series  2 and  3 tests  conducted  under 
this  program,  two  titanium  samples  ignited  and  were  consumed. 
The  conditions  for  these  tests  v/ere  1500°K,  1 atm  ©2  and 
1700°K,  .5  atm  02*  Examination  of  all  available  data  on 
sample  response  and  facility  performance  parameters  revealed 
that  in  each  case,  the  induction  heater  plate  current,  a 
measure  of  power  delivered,  increased  spontaneously  at  the 
time  the  sample  ignited.  It  cannot  be  stated  with  certainty 
if  the  ignition  process  proceeded  and  caused  the  abnormal 
induction  heater  response  or  if  the  induction  heater 
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FIGURE  2-4.  OXIDATION  DATA  TRENDS  AT  1500 °K 


malfunctioned  and  caused  the  sample  to  ignite.  Such  heater 
malfunction  is  not  likely  but  possible.  Repeats  of  the  test 
conditions  where  ignition  was  observed,  and  in  fact,  exceeding 
these  conditions  deliberately,  did  not  result  in  ignition 
supporting  the  heater  malfunction  explanation. 

2 . 3 Particle  Impact  Test  Program 

The  Phase  III  impact  tests  were  conducted  at  the  SAI 
Santa  Ana,  California  laboratory  using  the  two-stage  gas 
gun.  The  facility  is  shown  schematically  in  Figure  2-5. 

The  titanium  disc  specimen  with  a Type  E thermocouple  welded 
to  the  back  {Figure  2-2)  is  held  in  a sample  holder  and  the 
induction  coil  is  wrapped  around  the  holder.  A gas  supply  line 
(either  O2  or  H2O  vapor)  is  oriented  to  impinge  directly  on 
the  test  sample  providing  an  oxidizing  environment. 

The  test  procedure  consisted  of:  1)  initiating  gas 
flow  and  establishing  a 1/2  atmosphere  environment  of 
the  desired  gas;  2)  initiating  induction  heater  and  adjusting 
power  supply  to  attain  the  prescribed  temperature  (1200°K) ; 

3)  initiate  gun;  and,  4)  shut-down  heater  and  gas  flow.  The 
temperature  history  recorded  using  a thermocouple  was  used 
to  evaluate  if  catastrophic  ignition  occurred  after  the 
particle  impact. 

The  temperature  history  of  sample  was  monitored  during 
heat-up  and  impact.  For  all  test  conditions,  that  is,  inert 
and  oxidizing  atmospheres  and  single  and  multiple  glass  and 
single  ice  impacts  into  titanium  samples  at  1200°K,  no 
measurable  temperature  damage  resulted.  The  matrix  of  the 
Phase  III  tests  is  shown  in  Table  2-2.  The  impacts  produced 
craters  on  the  titanium  surface  and  spalled  off  the  oxide  layers 
over  large  areas.  As  expected,  the  heavier  glass  particles 
produced  more  damage,  that  is  larger  craters  (with  glass  residue 
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in  the  crater  bottom)  and  larger  spall  regions  than  the  ice 
particle.  A comparison  of  two  typical  craters  resulting 
from  ice  and  glass  impact  is  shown  in  the  post  test  photograph 
of  Figure  2.6. 
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3.0  AEDC  DET  EROSION  AND  OXIDATION  TEST  SERIES 

3.1  Background 

The  test  results  described  in  the  preceding  section 
show  the  effect  of  oxidizing  environments  (H2O  vapor  and  O2) 
on  titanium  under  both  non-erosive  and  single  impact  test 
conditions  using  both  non-reactive  (glass)  and  reactive 
(ice)  particles.  This  data  is  fundamental  to  the  develop- 
ment of  a predictive  model  for  titanium  response  as  will  be 
discussed  in  Section  4.0.  However,  the  test  conditions  are 
not  representative  of  the  multiple  impact,  hypersonic  flow 
environment  expected  for  a missile  structure  in-flight. 
Consequently,  a test  series  was  planned  which  would  use  the 
AEDC  Dust  Erosion  Tunnel  (DET)  to  evaluate  titanium  response 
to  multiple  particle  impacts  in  a hypersonic  flow  environ- 
ment. A facility  modification  was  implemented  to  introduce 
liquid  water  simultaneously  with  dust  (magnesium  oxide) 
particles  in  an  attempt  to  simulate  the  nuclear  cloud  water 
and  dust  environment. 

Extensive  tests  of  titanium  nosetip  models  have  been 
conducted  in  the  DET  in  previous  programs  (References  1-5) . 
These  tests  have  been  conducted  to  evaluate  increased  heat 
transfer  in  particle  environments.  Important  mechanisms 
include  particle  kinetic  energy  transfer,  particle-f lowf ield 
interaction  convective  heat  transfer,  and  high  temperature 
chemistry  heating  augmentation  due  to  surface  oxidation. 
These  effects  have  been  adequately  modeled  for  dust  part- 
icle/air environments.  However,  the  effect  of  water,  in  the 
form  of  ice,  liquid,  or  vapor  in  conjunction  with  a dust/air 
environment,  on  titanium  has  not  been  addressed.  The  con- 
cern was  that  water  could  increase  the  surface  chemistry 
effect  on  titanium  under  erosive  conditions  and  should  be 
considered  in  the  design  of  titanium  missile  structures. 

This  AEDC  DET  test  program  was  structured  to  evaluate  this 
concern . 
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Test  Objectives 

A test  program  was  planned  for  the  AEDC  DET  to  in- 
vestigate the  effects  of  water/air  mixtures  on  the  oxidation 
response  of  titanium  under  multiple  dust  particle  condi- 
tions. Although  the  DET  is  used  routinely  for  dust  erosion 
testing,  water  had  not  been  injected  prior  to  this  program. 

A facility  modification  was  implemented  to  inject  liquid 
water  either  separately  or  in  combination  with  dust  (MgO) 
particles.  The  goal  was  to  achieve  liquid  water  droplet 
impacts  on  the  hemispherical  titanium  models  which  had  been 
heated  to  high  temperature  by  the  hypersonic  air  and  dust 
flow.  This  test  program  was  the  first  such  utilization  of 
the  AEDC  DET  and  was  necessarily  exploratory  in  nature. 

3.3  Test  Program  Description 

3.3.1  Facility  Description 

The  tests  were  conducted  in  the  Arnold  Engineering 
Development  Center,  Propulsion  Wind  Tunnel  Facility,  Dust 
Erosion  Tunnel  (DET) . A description  of  the  DET  is  given  in 
Reference  6.  All  test  runs  were  made  at  the  127.5-in. 
nozzle  station,  where  the  exit  diameter  is  8.3  in.  and  the 
nominal  Mach  number  is  M = 6.6  to  7.5,  depending  on  oper- 
ating conditions.  A facility  modification  was  designed  and 
implemented  to  inject  liquid  water  simultaneously  with  the 
dust  particles.  A schematic  of  the  water/dust  injection 
strut  is  shown  in  Figure  3-1.  A separate  water  injection 
system  was  designed.  The  water  injection  line  is  shown  in 
Figure  3-1  next  to  the  dust  injection  line.  Various  ori- 
fices on  the  water  injection  line  were  evaluated  to  obtain 
the  desired  flow  conditions. 

3.3.2  Water  Orifice  Calibration 

A series  of  bench  tests  were  performed  to  calibrate 
the  water  orifice.  Four  designs  for  the  orifice  were  eval- 
uated. The  objective  of  the  calibration  tests  was  to  obtain 


a water  jet  with  minimum  divergence.  This  minimizes  the 
tendency  of  the  injected  water  to  atomize  in  the  throat 
region  of  the  AEDC  DET  and  maximizes  the  possibility  of 
water  impinging  on  the  model  at  the  test  section.  Test  data 
includes  photographic  records  of  the  water  jet  exiting  the 

orifice  and  a calibration  curve  of  injected  water  flow  rate 

• 

My  Q (gal/min)  versus  pressure  differential  AP  (psia) . The 
test  procedure  was  as  follows: 

1)  The  water  hopper  was  charged  with  a weighed 
amount  of  water  (distilled) . 

2)  The  hopper  was  pressurized  with  N2  ^rid  the 
orifice  flow  recorded  on  film  for  nominal 
AP's  of  200,  500,  and  800  psi. 

3)  After  each  flow  period  (1  minute  for  0.062 
orifice,  and  2 minutes  for  the  rest)  at 

each  AP,  the  water  was  carefully  drained  from 
the  system  and  weighed. 

4)  The  hopper  was  re-charged  and  the  sequence 
was  repeated  for  12  data  points. 

Table  3-1  summarizes  the  12  bench  tests  indicating  AP, 

« 

orifice  design  and  measured  My  q.  Figure  3-2  is  a schematic 
of  the  orifice  design  for  scenes  4-9.  The  0.02  inch  count- 
erbored  design  was  ultimately  chosen  for  testing  in  Runs  5- 
8.  Scenes  1-3  and  10-12  correspond  to  an  early  straight- 
through  design  shown  in  Figure  3-3  and  used  in  Runs  1-4, 
where,  as  indicated  in  Figures  3-4  - 3-7,  significant  diver- 
gence of  the  water  jet  was  obtained.  This  may  explain  the 
reason  for  no  apparent  water  impinging  on  models  in  early 
runs  (1-3)  for  low  injection  rates  (see  discussion  in  Sec- 
tion 3.4). 
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SCENE  AP,  PSI  DIA,  in  BORED  EXIT  INLET  RATE,  GPM 

1.34 
2.07 
2.68 
0.56 
0.88 
.1.18 
0.16 
0.25 
0.33 
0.15 
0.22 
0.28 


FIGURE  3-2.  CONICAL  INLET,  COUNTERBORED  EXIT  DESIGN 


FIGURE  3-3.  STRAIGHT  THROUGH  DESIGN 


Scene  2 
AP  = 500  psi 

~ 2.07  gal/min 


AP  = 200  psi 
Mh20  " 1.34  gal/min 


Scene  3 
AP  = 800  psi 
Mjj  Q = 2.68  gal/min 


FIGURE  3-4.  WATER  CALIBRATION  FOR  0.062  IN.  ORIFICE  (SCENES  1-3 
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0.56  gal/min 


0.88  gal/min 


1.18  gal/min 


j Figure  3-8  is  the  final  calibration  curve  for  each 

design.  Also  plotted  are  actual  values  for  test  runs.  Dif- 
ferences apparently  result  due  to  a difference  in  exit 
j orifice  pressure  from  static  throat  pressure  under  dynamic 

flow  conditions. 

3.3.3  Models 

The  test  models  were  6A1-4V  titanium  consisting  of 
eighteen  2.0  inch  diameter  (see  Figure  3-9  for  pre-test 
photos)  and  four  3.0  inch  diameter  models  (see  Figure  3-10 
for  pre-test  photos) . The  hemispherical  models  all  had  a 
0.25  inch  skirt  as  shown  in  Figure  3-11  for  attachment  to 
the  model  holder  shown  in  Figure  3-12.  Each  hemisphere  was 
instrumented  with  eight  platinum-platinum,  10  per  cent 
rhodium  thermocouples  welded  to  the  backside  at  the  loca- 
tions shown  in  Figure  3-11.  Table  3-2  presents  the  pre-test 
skin  thicknesses  for  each  model  and  thermocouple  location. 

3.3.4  Instrumentation 

In  addition  to  the  normal  DET  facility  measurements, 
primary  measurements  on  this  program  consisted  of  test 
section  total  pressure  and  total  temperature,  backface  model 
temperatures,  dust  and  water  injection  rates,  water  content 
in  the  test  gas,  and  motion  picture  coverage. 

A pressure  probe  was  used  to  measure  the  pitot  or 
normal  shock  recovery  pressure.  This  probe  was  a hemis- 
phere-cylinder with  one  inch  diameter  (Figure  3-13a)  and  a 
1/16  inch  diameter  orifice  at  the  sphere  stagnation  point. 
The  total  temperature  probe,  shown  in  Figure  3-13b,  used  a 
sheathed  Pt-Pt,  13  per  cent  Sh  thermocouple  with  a vapor- 
deposited  junction  shielded  by  a cylindrical  tube  made  of 
magnesia-stabilized  zirconia.  The  pressure  probe  and  temp- 
erature probe  outputs  were  recorded  on  conventional  strip 
chart  recorders.  Model  backside  thermocouple  outputs  were 
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FIGURE  3-9a.  PRE-TEST  PHOTOS  OF  2 . 0 INCH 
MODELS 


FIGURE  3-9f.  PRE-TEST  PHOTOS  OF  2.0  INCH  MODELS 
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FIGURE  3-lOa.  PRE-TEST  PHOTOS  OF  3.0  INCH  MODELS 


FIGURE  3-lOc.  PRE-TEST  PHOTOS  OF  3.0  INCH  MODELS 


FIGURE  3-11.  HEMISPHERE  MODEL  DESCRIPTION 
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OF  MODEL  HO 
POT  WELDED 


1/16  OD  Orifice  Tube 


a.  Pressure  Probe 


-Type  R Thermocouple  (38  gauge  wire)  in  0.063  in.  diam 
.Sheath.  Sputtered  Thin  Film  Junction 

\ ^4  Holes  90  deg  apart,  0.063  in.  diam. 

\ \ ^ No.  78  Sauereisen  Cement 

\ \ \ ^ Magnesia-Stabil ized  Zirconia, 

\ \ \ \ 1/4  00  X 3/16  ID 


f—  0.47  — sj  I 

0.75  

|« 0.97  »| 
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Not  to  Scale  - All  dimensions  in  inches. 


b.  Temperature  Probe 


FIGURE  3-13A&B.  PRESSURE  AND  TEMPERATURE  PROBES 
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recorded  on  magnetic  tape  using  an  analog-to-digital  con- 
verter, reduced  to  engineering  units,  then  input  to  an  IBM 
370/165  computer  for  execution  of  heat  transfer  calcula- 
tions. Initially,  attempts  were  made  to  measure  the  in- 
jected water  flow  rates  with  both  a turbine-type  flowmeter 
and  by  measured  pressure  differentials  across  the  orifice 
converted  to  volume  flow  by  a calibration.  These  efforts 
were  not  successful  because  the  turbine  meter  was  excep- 
tionally susceptible  in  electrical  noise  pickup  from  the  arc 
heater  and  the  AP  calibrations  could  not  be  performed  under 
the  same  conditions  as  the  orifices  were  subjected  to  during 
a test  run.  Beginning  with  Run  3,  the  water  was  weighed 
before  and  after  the  test  run  and  the  flow  rate  was  assumed 
constant  over  the  total  water  injection  time.  Dust  particle 
mass  flow  rates  were  also  computed  by  weighing  of  the  dust 
initially  put  into  the  dust  hopper  and  the  amount  remaining 
after  test  completion,  then  averaging  this  mass  difference 
over  the  total  dust  flow  time.  Three  motion  picture  cameras 
were  used,  one  for  overall  test  surveillance  and  two  for 
closeup  coverage  of  the  test  models.  These  cameras  were 
operated  at  50,  200,  and  400  fps.  A television  system  was 
used  to  observe  model  response  during  the  test  and  the 
viewing  screen  was  located  in  the  control  room. 

A gas  sampling  system  was  used  during  two  runs  to 
determine  the  water  content  of  the  test  gas.  The  system 
consisted  of  a 2 inch  diameter  stainless  steel  hemisphere 
with  a 1/16  inch  diameter  sampling  orifice  at  the  stagnation 
point,  two  stainless  steel  bottles  of  300  cc  capacity  each, 
and  several  valves  and  connecting  stainless  steel  tubing 
(see  Figure  3-14) . The  bottles  were  evacuated  and  dried 
prior  to  the  tests.  During  operation,  the  sampling  line 
from  the  model  was  purged,  then  the  remotely  operated  valves 
were  opened  for  about  10  seconds  each.  The  remotely  actuated 
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FIGURE  3-14.  GAS  SAMPLER  MODEL  (CONTINUED) 


valves  were  then  closed,  and  after  the  test  run,  the  hand 
valves  were  closed  and  the  bottles  sent  to  the  AEDC  Chemical 
Laboratory  for  analysis. 

3.3.5  Test  Procedures  and  Conditions 

The  sequence  of  events  for  a typical  test  run  was  as 
follows;  Prerun  calibrations  were  performed.  Specified 
tunnel  conditions  were  established  by  energizing  the  arc- 
heated  tunnel  driver  followed  by  20  seconds  for  stabili- 
zation of  cooling  water  temperatures.  With  the  tunnel  on- 
line and  ready,  all  data  recording  equipment  and  all  motion 
picture  cameras  were  started.  The  diagnostic  probes  and  the 
test  models  were  then  injected  into  the  test  stream  in 
specific  sequence  with  pre-programmed  periods  of  aerodynamic 
heating  and/or  dust  and  water  exposure.  After  all  models 
had  been  cycled  through  the  flow,  the  tunnel  driver  was  de- 
energized; then  the  postrun  calibrations  were  performed. 

Table  3-3  summarizes  the  AEDC  DET  test  conditions  for 
the  eight  runs  in  this  test  program  . Dust  particle  velo- 
cities and  densities  were  based  on  previous  AEDC  test  cali- 
bration data. 

3.4  Test  Results 

3.4.1  Test  Philosophy 

The  objective  of  the  DET  tests  was  to  demonstrate  the 
effects  of  water  on  titanium  surface  chemistry  under  particle 
erosion  conditions.  The  general  approach  to  accomplish  this 
objective  was  to  perform  tests  on  titanium  models  under 
varying  conditions  of  clear  air,  dust,  water  and  dust/water 
environments.  Comparisons  of  the  performance  of  these  model 
responses  were  made  to  identify  relative  effects. 

In  each  run,  9 stings  were  available  for  models  or 
instrumentation.  Typically,  these  stings  were  occupied  by 
one  pressure  probe,  one  temperature  probe,  four  models,  and 
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witness  bars  or  empty  stings.  For  the  high  enthalpy  runs, 
no  temperature  probe  was  used  since  it  would  not  survive. 

In  several  runs,  the  gas  sampler  was  used  to  measure  water 
content  in  the  test  gas. 

Primary  data  was  the  model  backface  temperature 
response.  Attention  was  focused  on  the  stagnation  region 
response.  Comparisons  of  this  temperature  data  and  reduced 
heat  transfer  data  for  the  stagnation  point  were  used  to 
determine  the  effect  of  water  on  the  titanium  response. 

This  test  program  encountered  a number  of  problems  in 
the  initial  stages  of  testing  related  to  facility  cali- 
bration and  operation  which  were  ultimately  corrected.  Part 
of  the  difficulties  were  associated  with  the  injection  of 
water  which  had  never  been  attempted  previously.  It  was 
especially  difficult  to  positively  confirm  the  presence  and 
state  of  the  water  at  the  test  section.  Witness  bars, 
motion  picture  coverage,  and  ultimately  the  gas  sampler 
instrument  were  used  to  confirm  the  presence  and/or  amounts 
of  water  at  the  test  section.  Originally,  four  runs  were 
planned  with  a required  expansion  to  eight  runs  to  accomp- 
lish the  objectives  of  the  program. 

The  detailed  test  data  results  for  Runs  2-8  are 
presented  in  Appendix  A.  Run  1 was  for  calibration  purposes 
and  is  not  discussed  in  this  report.  The  basic  data  shown 
in  Appendix  A includes:  (1)  post-test  model  photographs; 

(2)  total  temperature  probe  trace  (if  available);  (3)  model 
backface  temperature  measurements  for  all  thermocouple  loca- 
tions; and,  (4)  reduced  heat  transfer  data  for  all  thermo- 
couple locations.  A discussion  of  the  stagnation  point 
results  of  Runs  2-8  is  presented  in  the  following  section. 

3.4.2  Discussion  of  AEDC  PET  Runs 

Figures  3-15  - 3-21  present  the  results  of  AEDC  DET 
Runs  2-8.  The  backface  stagnation  point  thermocouple  temper- 
ature-time response  and  reduced  heat  transfer  are  shown  for 


each  of  the  models  tested.  Runs  2-6  models  were  2-inch 
diameter  and  Runs  7-8  models  were  3-inch  diameter  titanium 
{6A1-4V)  hemispheres  of  varying  ski  i thicknesses  --  model 
numbers  and  stagnation  point  skin  chicknesses  are  noted  on 
the  figures.  Test  conditions  including  the  dust  and  water 
environments  are  shown.  In  most  cases,  stagnation  pressure 
and  temperature  were  measured  with  probes  except  for  Run  2, 
where  the  temperature  probe  did  not  operate  correctly  and 
the  reported  value  is  based  on  an  energy  balance  and  Runs  7 
and  8 which  were  high  enthalpy  where  the  probe  would  not 
survive . 

Figures  3-15  and  3-16  show  the  results  for  the  ini- 
tial Run  2 with  nominal  flow  conditions  of  1000  psia  total 
pressure  and  800  Btu/lb  total  enthalpy  and  a low  water 
injection  rate.  These  results  show  the  expected  clear  air 
and  dust  response  including  the  high  temperature  ( ■2200°R) 
chemistry  effect  based  on  previous  studies  (Reference  4) . 

No  increased  effect  of  water  is  shown.  The  higher  temper- 
atures for  the  dust/water  model  Ti-6  compared  to  the  dust 
model  Ti-1  result  from  a thinner  skin  only.  In  fact,  no 
effect  of  water,  even  a cooling  effect,  is  shown.  Figure  3- 
16  shows  the  reduced  heat  transfer  rate  at  the  stagnation 
point.  This  reduced  data  is  characterized  by:  (1)  an 
initial  transient  increase  to  the  '\^1000-1500°R  regime;  (2)  a 
decreasing  heat  transfer  with  increasing  temperature  char- 
acteristic of  a convective  heating  trend  in  the  'Vl500-2000°R 
regime;  (3)  an  increasing  heat  transfer  with  increasing 
temperature  in  the  >2200°R  regime  characteristic  of  a dom- 
inating surface  chemistry  effect;  and,  (4)  periodic  eratic 
behavior  due  to  electrical  noise  in  the  thermocouple  output. 
Typically,  2000‘’R  is  an  excellent  temperature  to  compare 
relative  heat  transfer.  Figure  3-16  shows  no  increased 
effect  of  water  and,  in  fact,  possibly  a reduced  heat  trans- 
fer for  the  dust  and  water  model.  This  run  did  not  have  the 
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MODEL  BACKFACE  TEMPERATURE 


MODETi  TI-5  CLEAR  MODEL  TI-7  WATER 

Time  = 16.06  to  22.60  Sec  Time  = 42.97  to 
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gas  sampler  in  operation,  so  no  direct  evidence  exists  to 
prove  that  water  actually  impinged  on  the  model.  This  run 
was  the  original  0.020  inch  straight-through  water  orifice 
design  which  was  shown  in  Section  3.3.2  to  atomize  water  on 
injection  in  calibration  tests;  and  thus,  no  water  may  have 
impinged  on  the  model! 

The  conclusion  at  this  point  was  that  a large  un- 
certainty existed  as  to  the  presence  of  water  at  the  test 
section.  The  recommendation  was  to  reduce  the  total  en- 
thalpy (and  thus,  the  tendency  for  water  vaporization)  and 
increase  the  mass  flow  of  H2O  to  attempt  to  retain  the  water 
in  the  liquid  state  at  the  test  section. 

Figure  3-17  presents  the  results  for  the  Run  3 nom- 
inal 1000  psia,  500  Btu/lb  conditions.  This  run  again  used 
the  original  0.020  inch  straight-through  orifice  design  but 
more  water  was  injected  than  on  Run  2.  In  this  case,  there 
is  some  evidence  of  a slight  cooling  effect  of  the  water. 

The  dust  only  model  shows  a chemistry  effect  above  2000°F  as 
expected.  No  chemistry  effect  due  to  the  water  is  evident 
up  to  2000°F  in  this  run.  The  conclusion  at  this  point  was 
again  a lack  of  firm  evidence  that  water  was  present  in  the 
test  section.  The  recommendation  was  to  increase  the  size 
of  the  orifice  to  obtain  higher  water  mass  flow  rates  and  to 
reduce  the  total  pressure  to  maximize  the  potential  survival 
of  the  water  during  acceleration  by  the  air  flow. 

Figure  3-18  presents  the  results  for  the  Run  4 nom- 
inal 500  psia,  500  Btu/lb  conditions.  The  0.0625  inch 
straight-through  orifice  design  was  used  to  inject  a larger 
quantity  of  water.  The  clear  air  and  dust  only  response  was 
as  expected.  However,  the  water  only  and  dust/water  models 
experienced  massive  quenching.  These  models  did  not  reach 
temperatures  where  chemistry  would  be  expected  to  be  a 
factor . 
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FIGURE  3-17.  STAGNATION  POINT  RESULTS  FOR  AEDC  DET  RUN 
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FOR  AEDC  DET  RUN 


The  reconunendation  at  this  point  was  to  perform  the 
extensive  bench  calibration  tests  of  the  water  orifice  as 
described  in  Section  3.3.2  to  arrive  at  a better  orifice 
design  prior  to  repeating  the  previous  test  conditions.  In 
addition,  the  gas  sampler  was  designed  and  implemented  to 
quantify  the  water  content  of  the  gas. 

Figu  ‘es  3-19  and  3-20  present  the  results  for  the  Run 
5 nominal  500  psia,  500  Btu/lb  conditions.  This  run  used 
the  0.020  inch  counterbored  orifice  design  to  attain  a lower 
water  injection  rate  and  to  reduce  the  dispersion  of  the 
water  jet  as  it  emerged  from  the  orifice.  This  run  used  the 
gas  sampler  model.  The  primary  purpose  for  collecting  gas 
samples  was  to  measure  the  total  water  content  in  the  sample 
to  allow  a comparison  of  the  water/air  mass  ratio  of  the 
sample  to  the  water/air  ratio  determined  by  flow  rate  mea- 
surements. Two  samples  were  collected  during  each  run:  one 
sample  was  analyzed  for  total  water  content  and  the  other 
for  species  components.  Water  content  was  determined  by 
Karl-Fisher  titration  and  component  analysis  by  use  of  a gas 
chromatograph.  The  results  of  the  water  content  and  com- 
ponent analyses  are  given  in  Table  3-4.  For  Run  5,  the 

-2 

water/air  mass  ratio  of  injected  water  was  1.54X10  , and 

the  water/air  ratio  of  the  sample  was  6.54X10~^  or  42.5 
percent  of  the  injected  value.  The  test  procedure  was  also 
changed  for  Run  5 compared  to  preceding  runs.  Here,  the 
models  with  water  were  preheated  with  clear  air  or  dust  to 
high  temperature  prior  to  water  injection. 

Figure  3-20  shows  the  heat  transfer  rate  at  1500°R  to 

2 2 
be  reduced  from  20  Btu/ft  sec  for  clear  air  to  5 Btu/ft 

sec  for  the  water  only  model  (water  was  introduced  at  '^^1300°R)  . 

Note  that  dust  was  inadvertently  turned  on  and  accounts  for 

the  high  temperature  "turn-around"  in  heat  transfer  rate  for 
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the  clear  air  model.  Examination  of  the  dust  and  the  dust 
and  water  model  responses  also  shows  a cooling  effect. 

Water  was  introduced  at  the  end  of  the  dust  only  model  run 
and  accounts  for  the  significant  heat  transfer  reduction  at 
'v2200°R.  Water  was  introduced  at  'v2000°R  on  the  dust  and 
water  model  and  accounts  for  the  significant  heat  transfer 
reduction  observed.  No  chemistry  effect  of  water  is  ob- 
served for  these  models  which  experienced  bulk  temperatures 
to  'v2000‘’R. 

The  preliminary  conclusion  at  this  point  was  that, 
since  water  was  confirmed  at  the  test  section,  there  is  no 
net  effect  of  water  on  titanium  surface  chemistry,  at  least 
up  to  'V'2000°R  temperatures.  The  recommendation  was  to 
simulate  higher  total  enthalpies  to  confirm  or  deny  this 
preliminary  conclusion. 

Figures  3-21  and  3-22  present  the  results  for  the  Run 
6 nominal  1000  psia,  800  Btu/lb  conditions.  The  0.020  inch 
counterbored  orifice  design  was  used.  This  run  was  a repeat 
of  Run  2 except  the  dust  concentration  is  lower  and  the 
water  injection  rate  is  higher.  This  run  also  had  the  gas 

sampler  which  showed  that  water  was  impinging  on  the  model. 

_2 

For  Run  6,  the  water/air  ratio  injected  was  1.99X10  , and 

_2 

the  water/air  ratio  for  the  sample  was  2.06X10  or  103.4 
percent  of  the  injected  value.  These  results  for  the  5 
degree  off-stagnation  thermocouple  location  show  the  ex- 
pected surface  chemistry  effect  for  the  dust  only  model 
above  v2200°R.  The  effect  of  adding  water  is  shown  to  be  a 
net  cooling  effect  for  both  the  water  only  and  dust/water 
model.  The  preliminary  conclusion  was  that  water  does  not 
augment  the  surface  chemistry  for  this  intermediate  total 
enthalpy  flow  condition.  The  recommendation  was  to  conduct 
a high  enthalpy  flow  run  to  further  confirm  or  deny  this 
preliminary  conclusion. 
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MODEL  BACKFACC  TEMPERATURE 


MODEL  TI-12  CLEAR  MODEL  TI-2  WATER 

Time  = 29.86  to  34.30  Sec  Time  = 49.55  to  53.40  Sec 


k 

I 

I 

1 

Figures  3-23  and  3-24  present  the  stagnation  point 
backface  temperature  and  surface  heat  transfer  results  for 
the  Run  7 nominal  1000  psia,  1800  Btu/lb  conditions.  The 
.020  inch  counterbored  orifice  design  was  used.  The  gas 
sampler  could  not  be  used  to  confirm  the  presence  of  water 
since  it  would  not  survive  the  high  temperature  conditions. 
The  water  only  model  indicates  a net  cooling  effect  when 
water  is  initiated  at  '^'2500°R  compared  to  the  clear  air 
model  response.  The  dust  only  model  illustrates  the  ex- 
pected surface  chemistry  effects  above  'v2200°R.  Unfor- 
tunately, on  the  dust/water  model,  the  water  was  injected 
too  late  ('''3000‘’R)  to  obtain  extensive  information.  How- 
ever, a cooling  effect  is  suggested  at  these  high  temper- 
atures. The  recommendation  was  to  repeat  this  run. 

Figures  3-25  and  3-26  present  the  Run  8 nominal  1000 
psia,  1800  Btu/lb  results  for  the  repeat  of  Run  7.  Test 
conditions  were  essentially  identical  to  Run  7 except  that 
2.0  inch  models  were  used.  Again,  the  gas  sampler  could  not 
be  used.  The  results  show  the  same  increased  surface  chem- 
istry from  dust  above  '\^2200°R.  A decreased  heat  transfer 
I effect  of  water  is  shown  in  Figure  3-26  for  both  the  water 

j and  dust/water  models  compared  to  the  clear  air  and  dust 

j models,  respectively.  No  increased  effect  of  water  is  shown 

j for  bulk  temperatures  up  to  '^'3000°R. 
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TITANIUM  RESPONSE  MODEL 


4 . 1 Formulation  of  the  Titanium  Oxidation  Model 

The  oxidation  of  titanium  has  been  the  subject  of  many 
experimental  studies.  Earlier  work  has  been  reviewed  by 
Hauffe^'^^  (1965)  and  Kofstad^^^  (1966);  more  recent  work  is 
briefly  discussed  by  Wolf ^ (1975) . At  pressures  of  the 
order  of  1 atm  and  temperatures  above  600-700 °C,  experiment 
has  shown  that  the  oxidation  is  initially  parabolic,  i.e., 
mass  gain  % (time)^.  The  parabolic  behavior  regime  is  rela- 
tively long;  the  duration  is  about  500  minutes  at  800 °C,  de- 
creasing to  8 minutes  at  1200  “C.  Thus  for  the  problem  at 
hand,  i.e.,  modeling  titanium  nosetip  response  to  nuclear 
cloud  particle  impact,  where  the  surface  obscuration  times 
are  of  the  order  0.1  to  1 sec,  the  oxidation  can  be  considered 
to  be  wholly  within  the  parabolic  regime.  During  the  parabolic 
regime,  two  oxygen  consuming  processes  are  observed  to  occur; 

(a)  dissolution  of  oxygen  in  the  metal,  and  (b)  formation  of 
an  oxide  scale.  Our  objective  here  is  to  review  the  litera- 
ture to  assemble  the  data  required  to  develop  an  engineering 
model  of  the  oxidation  process.  In  order  to  identify  the  data 
required,  it  is  helpful  to  present  an  introductory  discussion  of 
the  proposed  oxidation  model. 


4.1.1  Proposed  Oxidation  Model 

The  proposed  model  for  parabolic  oxidation  postulates 
simultaneous  formation  of  a compact  scale  and  dissolution  of 
oxygen  in  the  metal.  The  schematic  below  depicts  the  situa- 
tion. 
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This  model  was  first  proposed  by  Wagner^^^^,  and  has  been  used 
to  interpret  zirconium  oxidation  data  in  (12,13).  In  these  ap- 
plications the  model  as  been  used  in  a semi-empirical  fashion 
to  avoid  requiring  the  knowledge  of  the  oxygen  concentration 
at  the  outer  surface  of  the  scale,  Cg.  Instead  a parabolic 
growth  law  for  scale  is  postulated  in  the  form 

e = 2 Y 

where  the  constant  Y is  determined  from  experiment.  Follow- 
ing this  approach  the  next  step  is  to  solve  the  species  con- 
servation equation  in  region  2,  the  metal.  In  our  final  model 
for  engineering  use,  this  solution  will  be  obtained  numerically 
since  O2  “ T is  a function  of  both  x and  t.  However, 

for  the  purpose  of  evaluating  data  in  the  literature,  it  is 
sufficient  to  consider  isothermal  oxidation,  for  which  an 
analytical  solution  can  be  obtained. 
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Solution  for  Isothemal  Oxidation 


In  Region  2 we  must  solve 


9c  _ p 9^c 

at  2 


subject  to 

t = 0,  c = Cq,  assumed  constant  here 
X = e,  c = Cy,  independent  of  time 
X c = Cq  = 0 say 

For  this  moving  boundary  problem,  we  follow  the  method  of 
Neumann.  A solution  to  Equation  (2)  which  satisfies  Equations 
(3)  and  (5)  is 


c = A erfc  

2/UJi 


at  X = e,  using  Equation  (4), 


c = A erfc  

2 /Pjt 


A 


^u 

erfc  Y 


Thus 


c 


erfc 

erfc 


X 

2/V^ 

Y'^Ci/P2 


is  the  desired  concentration  profile. 


r 


f 

f-: 


Next  we  postulate  that  the  parabolic  growth  law  for  the 
scale,  Equation  (1)  is  associated  with  a linear  concentration 
profile  in  the  scale,  then 


0<x<e,  c=Ca--  (c^  - c^) 


(9) 


In  the  limit  of  V2  ■*  0 it  is  easy  to  see  that  Equation  (5)  is 
consistent  with  Equation  (1) , for  then 


P = -F  ('^a  - '^s) 


Pi 


Pi  , .. 


or  £ a /t 


The  total  oxygen  present  in  the  scale  and  metal  at  time  t is 


= / 


c dx 


00 

= J" c dx  + J c"d 


^ c dx  +y  cdx  with  e 
o'^  e 


2y/^ 


/ 


c dx  = ^ (c^  + c„)  e 


(10) 


00  C 

fcdx  = f 
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£ erfc  Y 


erfc  — ^ — I <jx 

2/V^  I 
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2/17^ 


J erfc  n dn 

e 


2 /t^2 1 


4p2t 


^2^  - ^1/^2 


UI|  - 


erfc 


erfc 


Y 


] 


Thus,  combining  Equations  (10)  and  (11), 


^02  = 2 '^s^  + 

- 2y 


u 


erfc  Y /Pj^/P2 


^2"  - Y P1/P2 

^ " © 


Next  we  eliminate  the  unknown  c^  in  favor  of  y;  if  we  assume 
this  is  particularly  easy,  then 


1 

I 


Cg  2 Y ^ 


and 


M02 


<•=3  ■ ‘^u)  2 Y + 

erfc  Y 


which  agrees  with  Equation  (21)  in  Reference  (5).  The  assumption 
Ca  = Cg  proves  to  be  sufficiently  good  to  make  Equation  (12) 
generally  valid. 

For  the  rate  of  oxygen  consumption  we  write 


dM02 

dt 


K,  t"^/^  + K t“^^^ 
f s 


(15) 


where  rate  constants  K^,  and  have  been  introduced  for 
the  total  oxidation,  oxide  film,  and  solution,  respectively. 

Thus 

Kf  = (Cg  - c^)  Y (16) 


K . c /?-2 

S U V IT  — — ^ 

erfc  Y ^Pi/P2 


The  preceding  analysis  shows  that  the  parameters  of  the  model 
are  the  concentrations  c,,  c_  and  c,,  and  the  diffusion  coef- 
ficients  P^^  and  p2*  1*^  order  to  use  the  model  for  engineering 

analysis  data  for  these  parameters  is  required. 

4.1.2  Survey  of  Available  Data 

4.1.2. 1 P2  The  Diffusion  Coefficient  of  Oxygen  in  Titanium 

Considerable  data  for  P2  for  unalloyed  titanium  exists, 
but  data  for  titanium  alloys  is  sparse.  The  most  useful  data 
is  summarized  below. 
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1.  Kofstad  et  al  (14)  derived  a value  for  V2  for  unalloyed 
titanium  neglecting  the  consequences  of  oxide  film  formation, 

^2  = 0.5  exp  ( - cmVs  (19) 

in  the  temperature  range  800  - 1200°c. 

2.  Shamblen  and  Redden  (15)  obtained  for  7i-6Al-2Sn-4Zr-2M 

o 

^2  " ^o  cm  Vs  (20) 

= 147  [O2  ppm]‘°''^^2 
= 062  cm^/s 

in  the  temperature  range  300  - 870°C.  These  authors  found 
that  O2  increases  somewhat  with  substitutional  alloying  con- 
tent . 

3.  Rosa  (9)  determined  P2  a-titanium  from  the  rate  of 

movement  of  the  a-6  interface  as 

P2(a)  = 0*778  exp  (-  i^-^)  cm  Vs  (21) 

for  commerically  pure  titanium  in  the  temperature  range  932  - 
1142®C.  Rosa  also  determined  ^2  for  S-titanium  from  micro- 
hardness measurements  as 


for  the  same  conditions.  Rosa  makes  favorable  comparisons  of 
his  results  with  earlier  data  in  the  literature,  (17,18,19). 

Figure  4-1  shows  a comparison  of  Equations  (18,  19,  20 
and  21) . Also  shown  is  a "corrected”  Kofstad  equation  reder- 
ived from  his  oxidation  rate  data,  as  follows: 


From  Equations  (13)  and  (17) , for  dissolution  only. 


t 


-1/2 


Mo2  = 2K3  t^/2 


or 


Mn  = ^ c P-t  = K t 
^2  IT  u 2 p 


For  c^  = 0.641  gm/cm  (derived  below) 


0.523  K Po  = K 
P 2 p 


0.523  P exp 
o 


51,000 

RT 


Ap  exp 


51,000 

RT 


Kofstad  obtained  Ap  = 3.2,  thus  Pq  = 3.2/5.23  = 6.12  cm^/s 
which  is  an  order  of  magnitude  higher  than  the  value  given  by 
Equation  (19).  We  conclude  that  there  is  a numerical  error 
in  Equation  (19) . We  also  note  that  we  would  expect  the 
Kofstad  value  of  P2  to  be  higher  than  those  of  Rosa  and  Shamblen 
since  in  the  former  case  it  was  assumed  that  all  the  oxygen 
consumed  went  into  solution,  i.e.,  the  Ti02  formation  was  ignored. 
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FIGURE  4-1.  DIFFUSION  COEFFICIENT  IN  TITANIUM  FROM  THE  LITERATURE 
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4. 1.2. 2 c^.  The  Concentration  of  Oxygen  in  Titanium  at  the 

Ti02-Ti  Interface 

If  oxygen  is  available  in  excess  then  c^  would  be  ex- 
pected to  be  saturation  concentration  of  O2  in  solution  in  Ti. 
Figure  4-2  shows  the  Ti-0  phase  diagram  where  it  is  seen  that 
aTi  is  saturated  at  0.143  weight  per  cent  oxygen,  essentially 
independent  of  temperature.  Rosa  (16) , in  analyzing  his  diffu- 
sion coefficient  data  used  a value  of  c^  = 0.665  gm/cm-^  which 
corresponds  to  0.148  weight  per  cent.  Shamblin  and  Redden  (15) 
used  a value  of  0.143  when  analyzing  their  diffusion  coeffic- 
ient data.  Kofstad  et  al.  (14),  based  on  hardness  measurements, 
conclude  that  the  oxygen  concentration  in  the  outer  layers  of 
metal  tends  towards  a limiting  value  corresponding  to  a com- 
position TiOQ  35,  i.e.,  0.10  weight  per  cent.  Considering 
the  difficulty  of  ma)cing  hardness  measurements  very  close  to 
the  surface,  we  view  this  result  as  being  not  inconsistent  with 
a true  interface  value  of  0.143.  However  we  also  note  that  in 
Reference  (20),  a saturation  oxygen  concentration  of  10.02  per 
cent  is  claimed  but  no  citation  or  basis  for  this  figure  is 
given . 

The  parameter  c^  does  play  a central  role  in  the  oxida- 
tion process.  If,  as  available  evidence  seems  to  indicate, 

Cy  readily  attains  the  equilbrium  saturation  value  for  oxygen 
pressures  in  the  range  of  interest,  then  the  rate  of  solution 
of  oxygen  will  be  independent  of  oxygen  partial  pressure, 
and  also,  will  not  be  affected  by  the  addition  of  water  vapor 
to  the  ambient  atmosphere.  Thus  any  affect  of  oxygen  partial 
pressure  on  the  total  oxidation  rate  must  be  an  effect  on  the 
rate  of  growth  of  the  oxide  scale. 
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TEMPERATURE 
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FIGURE  4-2.  PHASE  DIAGRAM  FOR  THE  TI-0  SYSTEM 
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TEMPERATURE 


4. 1.2. 3 1^1  f Cg,  Cg  and  y,  The  Parameters  Governing  the  Oxide 

Scale 

Three  parameters  are  required  to  characterize  the  rate 
of  growth  of  the  oxide  scale;  Cg  and  (c^  - c^)  , or 

Cg  and  Y»  The  concentration  Cg  is  simply  that  of  oxygen  in 
Ti02»  0.40  weight  per  cent.  However  no  data  for  (c^  - Cg) 
or  Y appear  in  the  literature.  In  fact,  although  the  Wagner 
model  previously  described  is  widely  ignored,  there  is  doubt 
as  to  whether  the  diffusion  across  the  oxide  film  can  be  des- 
cribed by  Picks  law  with  oxygen  concentration  as  the  driving 
potential.  Whether  the  oxide  scale  is  predominatly  an  electronic 
conductor,  and  whether  it  is  a p-  or  n-type  conductor,  depends 
on  temperature.  Kofstad  (9)  summarizes  existing  theory  for 
transport  across  metal  oxide  scales. 

For  modeling  purposes  we  require  only  a single  addi- 

1/2 

tional  parameter,  the  combination  f (^a  ~ ^s^  '<'^’1  • 

Thus  we  recommend  that  available  data  be  used  to  determine 
either  combination  without  explicit  considerations  of  the 
transport  processes  in  the  scale.  Two  possible  methods  are 
(1)  use  of  data  for  the  relative  amounts  of  oxygen  going  into 
solution  and  into  the  scale,  and  (2)  use  of  data  for  the  total 
oxidation  rate  constant  Kq  together  with  valuer  of  V2  and  c^/ 
and  Equation  (13-17).  The  former  method  is  particularly  simple 
but  unfortunately  contradictions  appear  in  the  literature. 

Kofstad  et  al.  (14)  report  a fraction  of  20%  of  oxygen  going 
into  the  scale  at  1173  K,  while  Stringer  (21)  reports  a value 
of  55%  at  1223  K.  Other  studies  have  not  presented  quantita- 
tive data,  but  suggest  that  most  of  the  oxygen  goes  into  solu- 
tion, particularly  at  higher  temperatures.  However,  even  if  we 
were  to  disregard  Stringer's  result,  further  data  over  a range 
of  temperature  would  be  required  to  determine  the  activation 
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energy  of  Application  of  the  latter  method  requires  selec- 

tion of  the  most  reliable  and  appropriate  data  for  and 
since,  again,  contradictions  do  appear  in  the  literature.  For 
example,  if  we  were  to  use  the  Kofstad  data  for  and  Rosa's 
data  for  V21  we  would  deduce  that  considerably  more  than  20% 
of  the  oxygen  was  going  to  form  scale,  in  variance  with  the 
value  reported  by  Kofstad. 

4. 1.2. 4 Data  for  Kc,  The  Total  Oxidation  Rate  Constant 


For  unalloyed  titanium  Stringer  (21)  compared  the  data 

of  his  own  and  several  previous  studies.  Figure  4-3  taken  from 

”2 

(21)  shows  his  result.  His  correlation  for  Kp  (=4K^)  is 

4K^2=  1.007  x 10“^  exp  (-  + 320  exp  (-  gmVcm^s  (23) 

and  concluded  that  V-^  has  an  activation  energy  of  2 5,4  00 
cal/mole  while  P2  has  an  activation  energy  of  61,800  cal/mole. 

It  is  of  interest  to  note  that  Kf/K^  calculated  from  Equation 
(23)  at  1223  K is  close  to  unity  and  not  55%  as  quoted  by 
Kofstad  (9)  as  being  obtained  by  Stringer. 

Kofstad  (9)  compared  the  data  obtained  by  Kofstad  et  al. 

(14)  with  previous  studies,  as  shown  in  Figure  4-4:  good  agree- 
ment is  apparent.  The  more  recent  data  is  that  of  Wolf  (lo). 

Figure  4-5  shows  Wolf's  data  for  "Type  1"  unalloyed  titanium 
with  Stringer's  correlation,  Equation  (22)  and  the  data  of 
Kofstad  et  al.(i4).  VJithin  experimental  error  there  is  excel- 
lent agreement  between  Wolf's  data  points  and  the  Stringer 
correlation.  Comparison  of  Figures  4-3,  4-4  and  4-5  leads 
to  the  conclusion  that  the  experimental  data  for  K^  is  reason- 
ably consistent  and  well  represented  by  Equation  (23). 
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. SUMMARY  OF  SPECIMENS  (TI6AL-4V)  USED  TO  COMPARE  PREDICTED 
REACTION  RATE  WITH  EXPERIMENT  (10) 


T. ’C 

1200  ICOO  900  800  700  600 


1-- 


6 7 8 9 10  U 12 

10*/T,  ’K 


Figure  4-4.  Total  Oxidation  Rate  Constant  [KOFSTAD] 


I 


i 

Pressure  Dependence 

Since  a large  number  of  investigations  for  unalloyed  i 

titanium  are  well  correlated  by  Equation  (23)  , it  is  possible 
j to  discern  if  Kq  is  dependent  on  oxygen  pressure.  Table  4-1 

below  lists  the  oxygen  pressures  considered.  We  conclude  that 
there  is  no  dependence  of  oxidation  rate  on  oxygen  pressure 
in  the  range  of  conditions  under  consideration.  Wolf  (10)  ob- 
tained one  data  point  at  400  torr  for  which  was  about  twice 

the  value  be  obtained  at  200  torr  and  the  same  temperature. 

He  thus  claims  that  "the  degree  of  reaction  is  very  nearly 
linearly  related  to  oxygen  pressure.  This  speaks  strongly 
against  any  type  of  solid-state  reaction  control..".  It 
seems  quite  unwarranted  for  Wolf  to  make  such  a claim  based  on 
one  data  point,  when  there  is  so  much  evidence  in  the  liter- 
ature to  the  contrary:  we  are  of  the  opinion  that  Wolf's  re- 
sult should  be  disregarded  at  this  time. 

Table  4-1.  Oxygen  Pressures  in  Investigations 
of  Unalloyed  Titanium 

INVESTIGATION 

Wolf  (10) 

Stringer  (21) 

Kofstad  et  al.  (14) 

Jenkins  (22) 

Davies  and  Birchenall  (23) 

Hurlen  (24) 


PRESSURES 

200  torr 
1 atm 
1 atm 
40-700  mm 
1 atm 
1-760  mm 
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Ti-6Al-4V  Alloy 

The  only  data  available  for  this  alloy  is  that  of 
Wolf  (10)  , which  is  also  shown  in  Figure  4-5.  The  oxidation 

rate  is  higher  than  that  for  unalloyed  titanium  which  is  in 
line  with  the  conclusion  of  Shambler  and  Redden  (15)  that  1*2 
increases  somewhat  with  substitutional  alloying  content. 

Given  the  scatter  in  the  experimental  data  there  is  no  basis 
for  concluding  that  the  activation  energy  differs  from  that  of 
unalloyed  titanium.  For  temperatures  above  1000  C Equation 
(23)  can  be  approximated 

4k ^ s 320  exp  (-  gm^/cm^s  (unalloyed)  (24) 

O K1 

2 

For  Ti-6Al-4V  the  values  of  measured  were  on  an  average 

60%  higher  than  Equation  (23)  , thus  we  recommend 

4Kq^  = 480  exp  (-  gm^/cm^s  (Ti-6Al-4V)  (25) 

4.1.3  Oxidation  of  Titanium  by  Steam 

At  the  temperatures  and  pressures  of  interest  steam  will 
be  only  slightly  dissociated.  Table  4-2  below  tabulates  equi- 
librium compositions  at  0.3  atm  pressure. 

Table  4-2.  Equilibrium  Compositions  of  Steam  at 
0 . 3 atm  ( 25 ) . 


COMPOSITION 

temperature”  k 

(MOLE  FRACTIONS) 

1200 

1400 

1500 

1700 

»2 

0 

.0001 

.0003 

.0007 

H2O 

1.0000 

. 9998 

.9995 

. 0089 

O2 

0 

0 

.0001 

.0003 

OH 

0 

0 

.0001 

.0002 
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When  titaniiom  reacts  with  steam  an  oxide  layer  forms  in  a very 
similar  manner  as  in  the  process  observed  with  oxygen;  thus,  it 
is  probable  that  the  H2O  diffuses  through  the  oxide  film  and 
dissociates  at  the  Ti02  - Ti  interface.  The  equilibrium  of 
such  a reaction  has  not  been  investigated  in  the  literature; 
also,  there  have  been  no  studies  of  H2O  diffusion  through  Ti02. 
The  process  may  be  further  complicated  by  disparities  in  the 
rate  of  diffusion  of  oxygen  and  hydrogen  into  titanium. 

Data  for  Kq 


Lohberg  and  Schleicher  (26)  have  investigated  the  titanium- 
steam  reaction  in  the  range  600-1000  C.  Figure  4-6  shows  their 
results  where  it  is  seen  that  a parabolic  steam  consumption 
curve  was  not  obtained,  instead 

»H20  “ 


a 


(27) 


The  activation  energy  for  steam  consumption  was  determined 
to  be  20,400  - 20,900  cal/mole.  In  (26)  the  authors  state 
that  Lutsch)cin  and  Iljin  (27),  in  an  earlier  study,  obtained 
an  activation  energy  of  19,100  cal/mole.  If  a parabolic  oxi- 
dation law  were  assumed  these  values  will  correspond  to  an 
activation  energy  for  V2  about  40,000  cal/mole,  somewhat 
lower  than  that  obtained  for  oxygen.  Figure  4-7  compares  the 
gas  consumption  at  1000 “C  with  that  for  oxygen,  as  given  by 
Equation  (24) . 
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Figure  4-6.  Titanium  Steam  Reaction 
Lohberg  and  Schleicher 
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Lainer,  Bai  and  Typsin  (28)  also  studied  the  titanium- 
steam  reactions  but  the  data  reported  is  sparse.  Oxidation 
rates  were  correlated  with  the  diffusion  coefficient  V-^  for 
oxygen  transport  in  a rutile  scale.  Values  for  V-^  formed  by 
oxidation  in  H2O  are  claimed  to  be  40-600  times  the  values 
obtained  for  oxidation  in  air.  No  rate  constants  for  the 
titanium-steam  reaction  are  reported. 

4.1.4  Heats  of  Oxidation  and  Solution  of  Oxygen  in  Titanium 

The  energy  evolved  per  gm  atom  of  reacting  oxygen  in 
titanium  for  different  end  products  has  been  calculated  by 
Mah  et.al.  (Bureau  of  Mines  Report  of  Investigation  5316) 
and  is  shown  in  Figure  4-8.  For  modeling  purposes,  it  is 
adequate  to  use  as  heats  of  oxidation  and  solution,  12024 
Btu/lb  O2  and  15,575  Btu/lb  O2  respectively. 
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Figure  4-8.  Energy  Produced  per  gm-atom  of  Reacting  Oxygen 
(16  GM  0)  IN  Ti  FOR  Different  End  Products 
and  Reaction  Temperature 


ISBEStaUALlTIJBACn^ 


4-22 


4.1.5  Hydrogen  Dissolution  in  Titanium  Review 


4. 1.5.1  Literature  Survey 

The  dissolution  of  hydrogen  in  titanium  differs  from 
oxygen  dissolution  in  that  it  is  a reversible  reaction,  and 
provided  the  temperature  is  high  enough  to  permit  finite 
rate  processes,  titanium-hydrogen  alloys  can  only  exist  in 
equilibrium  with  a gaseous  hydrogen  atmosphere.  Figure  4-9 
shows  the  equilibrium  concentrations  for  pure  titanium. 
Figure  4-10  shows  the  phase  diagram  and  Figure  4—11  is 
deduced  from  4-10. 


FIGURE  4-9.  EQUILIBRIUM  CONCENTRATIONS  FOR  PURE  TITANIUM 


Waisman,  Sines  and  Robinson  (29)  investigated  equilibrium 
for  Ti-6Al-4V  alloy  and  found  that,  at  a given  solid  phase 
hydrogen  concentration  and  temperature,  the  equilibrum  gas 
pressure  was  2.31  times  that  for  pure  titanium. 
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FIGURE  4-10.  PHASE  DIAGRAM  FOR  THE  HYDROGEN-TITANIUM  SYSTEM 
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FIGURE  4-11, 


TITANIUM  PHASES  AS  THEY  DEPEND  ON  TEMPERATURE  AND 
HYDROGEN  PRESSURE 
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Diffusion  coefficients  of  hydrogen  in  pure  titanium 
have  been  measured  by  Wasilewski  and  Kehl  (30)  as 


P.  = 1.95  X lO'^  exp  { (-6640  + 500)/RT}  cm^/s 

P 


= 1.8  X 10~^  exp  { (-12380+680)/Rt}  cm^/s 


which  are  considerably  larger  than  the  diffusion  coeffi- 
cients for  oxygen  in  titanium. 

The  heat  of  solution  of  hydrogen  in  titanium  has  been 
calculated  by  McQuillan,  and  is  reported  by  Skinner  et.al. 
(31)  as  -21,600  cal/mole  to  form  a solution  and  -27,800 
cal/mole  to  form  6 solution.  These  values  are  valid  for 
dilute  solutions  (less  than  10  atom  % H) . 


4. 1.5. 2 Sample  Numerical  Calculations 

Since  the  mechanism  and  equilibrium  of  the  titanium- 


steam  reaction  is  unknown,  it  is  only  possible  to  obtain  rough 
estimates  of  hydrogen  dissolution  in  titanium  exposed  to  steam. 
Since  the  diffusion  coefficient  for  hydrogen  in  titanium  is 
considerably  larger  than  that  for  oxygen,  a possible  approach 
is  to  simply  assume  that  the  hydrogen  goes  into  solution  in 


the  stoichiometric  proportion  indicated  by  the  H2O  formula. 


i.e.,  1 lb  H for  8 lbs  At  1200K  the  oxygen  reaction  rate 

is 


M 


02  = 


-1/2 


K = (4K^) 
O p 


1/2 


, , -30,900  , 

-22.4  exp  ( (1.987)  (1200) ' 


= 5.27  X 10“^  Ib/ft^s^^^ 


M = 5.27  X 10"^t"^/^  Ib/ft^s^/^ 
°2 


M 
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AH 


Heat  evolved 


M, 


(5.27  X 10 

0.83 


■5)  (^^P-^.OQP.)  (1.8)t'^^^ 

Btu/ft^s^/^ 


Then 


0.125  M 

0.659  X 10”^t~*^  Ib/ft^s^ 


Heat  evolved  = AH 

"2 

= (0.659  X io~^) (27,800)  (i 
= 0.165  t“^Btu/ft^s^ 


which  is  about  20%  of  the  heat  evolved  by  the  O2  reaction. 

However,  McQuillan  (32)  states  that  the  rate  of 
hydrogen  absorption  is  considerably  reduced  when  the  surface 
has  been  "contaminated"  by  oxygen.  No  quantitative  infor- 
mation on  the  effect  has  been  published,  but  Gulbransen  and 
Andrew  (33)  have  examined  the  corresponding  effect  in  Zir- 
conium in  some  detail.  A further  and,  perhaps,  more  im- 
portant consideration  is  that  the  bulk  of  the  water  vapor 
available  for  reaction  with  a titanium  shroud  results  from 
evaporation  of  water  droplets  or  ice  particles.  Let  us  now 
calculate  the  latent  heat  associated  with  a droplet  evapor- 
ation process  to  produce  the  hydrogen  dissolution  rate 
calculated  above.  The  amount  of  water  evaporated  to  sustain 
the  hydrogen  dissolution  rate  calculated  above  is 
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= (18/2)  (.659  X 

2 

= 5.93  X Ib/ft^s** 

and  since  the  latent  heat  is  approximately  1000  Btu/lb  the  rate 
at  which  heat  is  absorbed  is 

= (5.93  X 10"^)  (1000)  t"**  = .059t~‘*  Btu/ft^ 

v«p 

which  is  36%  of  the  heat  evolved  by  the  hydrogen  dissolution. 
Thus,  at  most  the  net  increase  in  heat  liberation  is  (1  - 0.36) 
(20)  = 12.8%  more  than  the  value  obtained  for  oxidation  in  air 
or  oxygen. 

4 . 1 . 5 . 3 Effect  of  Temperature  Level 

At  higher  temperatures  the  degree  of  dissociation  of 
H2O  will  be  greater  (see  Table  4-2) , and  the  diffusion  co- 
efficient of  H in  titanium  is  also  higher;  however,  the 
solubility  of  H decreases  rapidly  with  temperature  (see 
Figure  4-9)  and  the  net  effect  is,  for  example,  at  1500K  the 
rate  of  dissolution  of  H is  about  40%  less  than  at  1200K. 


= (81/2)Mjj 


4.1.6  Evaluation  of  Clemaon  Data  (3)  for  Oxygen  Consumption 
Rate;  Ti-6Al-4V  Alloy 

4. 1.6.1  Model  Parameters  Based  on  Weight  Change  Measurements 

Based  on  the  long  time  weight  change  measurements 
made  at  Clemson  (10)  we  proposed  Eq . (20)  to  give  the  rate 

constant  Kp  for  oxidation  of  Ti-6A1-4V  alloy,  viz. 

Kp  = 480  exp  ( — ) gm  /cm  s (24) 


Our  model  parameters  1^2'  ^ calculated  from  Kp 

as  follows.  Ignoring  moving  boundary  effects  and  postu- 
lating 80%  of  the  oxygen  goes  into  solution  while  20%  goes 
to  form  oxide,  we  can  derive 


2 2 
4k::  = K^ 
o p 


Ti02  2^,2  2 ^2  ^1 

■3713  V < * =u  4 «a-S> 


= 0 (26) 


^2  1/2  2 1/2 
2C  (^)  = (0.8)  (4K  ^) 

U iT  O 


f 


^2  1 /2  C +C  1 /T 

2C^  (-^)  t^/^  = 4(-a,2-^)  2 Y(P^t)^/2 


with  =40.2  Ib/ff 


C = 106.3  Ib/ft' 
s 


= 107.2  Ib/ff 

a 


P^io^  = 281  Ib/ff 


and  assuming  the  same  actiation  energy  for  and  ^2' 
equations  may  be  solved  to  give 
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Y = .0297 

= 2.03  ftVs 
p2Q  = 0.676  ft^/s 

4. 1. 6. 2 Comparison  with  Volumetric  O2  Consumption  Measurements 

Table  4-3  summarizes  the  pertinent  details  of  the  tests 
for  which  satisfactory  measurements  of  oxygen  consumption 
rate  were  made.  Turning  first  to  the  low  heating  rate  tests 
(0.5®C/s),  Figure  4-12  shows  a comparison  of  the  predictions 
of  our  computer  code  with  the  experimental  data  for  Mq  as  a 
function  of  time.  Since  the  prediction  lies  between  t^e 
experimental  results  for  the  two  types  of  Ti-6A1-4V  alloy, 
we  view  the  prediction  as  quite  satisfactory.  For  these  low 
heating  rates,  Wolf  et.al.  (34)  also  concluded  that  the 
parabolic  kinetics  implicit  in  the  model  are  satisfactory. 

Figure  4-13  shows  a comparison  at  a heating  rate  of 
8®C/s  where  the  comparison  is,  at  first  sight,  less  satis- 
factory. However,  a number  of  considerations  are  pertinent. 

(i)  In  the  Clemson  experiments,  there  were  apparent  diffi- 
culties in  recording  the  zero  time  point  accurately.  Table 
3 compares  the  time  that  should  have  elapsed  to  reach  1100‘’C 
(calculated)  with  the  actual  time;  the  discrepancies  are  40, 

45  and  30  seconds  for  the  three  specimens.  Thus,  in  Figure 
13  the  experimental  curves  have  been  adjusted  by  the  above 
time  displacements,  a rather  unsatisfactory  procedure.  (ii) 
Since  the  three  specimens  are  of  the  same  type  of  Ti-6Al-4V 
alloy,  the  differences  between  the  experimental  curves 
indicate  the  lack  of  precision  in  the  experiments  and  pos- 
sible variability  of  physical  phenomena  (e.g.,  cracking  of 
the  oxide  scale) . (iii)  The  data  points  indicated  on  the 
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FIGURE  4-13.  COMPARISON  OF  SAI  PREDICTIONS  AND  EXPERIMENTAL  RESULTS  FROM 
REFERENCE  10  (8°  c/s) 


experimental  curves  are  those  actually  recorded.  It  is 
clear  that  they  are  too  sparse  to  record  a peak  in  as 
sharp  as  indicated  by  theory.  Furthermore,  in  the  experi- 
ments there  were  problems  in  preventing  temperature  over- 
shoots at  the  end  of  the  heating  period.  Such  overshoots 
would  tend  to  broaden  the  observed  peaks  in  M_  . Figure  4-14 
shows  a comparison  at  a heating  rate  of  22°C/s;  again,  the 

sharpness  of  the  predicted  peak  in  is  in  evidence. 

^2 

In  our  review,  the  predictions  at  higher  heating 
rates  are  adequate  in  view  of  the  experimental  uncertainties. 
Thus,  we  do  not  support  the  view  of  Wolf  et.al.  (34)  that  at 
high  heating  rates  the  parabolic  kinetics  should  be  dis- 
carded in  favor  of  the  asymptotic  rate  law  of  Evans  (35). 
Furthermore,  we  are  also  doubtful  if  further  corrections 
should  be  made  as  suggested  in  (10)  , for  example,  the  sug- 
gested increase  in  effective  reaction  surface  area  during 
rapid  heating.  It  is  important  to  note  that  the  intro- 
duction of  an  ad  hoc  kinetics  expression  such  as  the  asymp- 
totic rate  law  would  preclude  rational  modeling  of  the 
oxidation  process : it  would  then  not  be  possible  to  account 
for  such  features  as  spatial  temperature  variations,  and 
mechanical  removal  of  the  oxide  film. 
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4.2  Analysis  of  Oxidation  Test  Data 

The  analytical  model  developed  in  Section  4.1  for 
oxidation  of  6AL-4V  titanium  in  oxygen  environments  was  used 
to  analyze  the  results  of  the  Phase  I SAI  oxidation  test 
data.  Analytical  predictions  for  the  weight  gain  as  a 
function  of  time  for  the  indicated  experimentally  measured 
temperature  history  are  shown  in  Figure  4-15  for  typical 
runs.  The  final  experimental  weight  gain  is  shown  for 
comparison . 

Important  conclusions  from  this  analysis  include: 

1)  weight  gain  during  the  transient  heat-up 
period  is  not  important 

2)  the  tests  were  essentially  isothermal 
temperature  oxidation  governed  by  the 
parabolic  law  as  developed  in  Section  4.1. 

The  total  experimental  data  base  is  shown  in  Figure  4-16  in 
terms  of  weight  gain  (AVJ)  normalized  by  square  root  of  time 
(reflecting  the  parabolic  law)  in  comparison  with  the  anal- 
ytical model  and  uncertainty  bands  based  on  information  in 
the  literature.  The  Phase  I SAI  data  falls  within  the 
uncertainty  bounds  of  the  model  and  generally  supports  the 
parabolic  oxidation  rate  law.  Within  the  context  of  this 
evaluation,  an  increased  oxidation  in  H2O  vapor  is  shown  up 
to  1700®K.  A pressure  effect  is  also  apparent,  especially 
at  1400‘’K,  although  the  1400°K  data  appear  anomalous  com- 
pared to  the  other  data  points. 

Application  of  the  oxidation  model  to  the  AEDC  DET 
data  shows  good  agreement  with  measured  heat  transfer. 

Figure  4-17  illustrates  the  Run  7 ti-22  model  clear  air 
response  stagnation  point  heat  transfer  data  in  comparison 
with  predictions.  The  dominant  convective  heating  trend 
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occurs  until  2200®R  where  surface  oxidation  heating  becomes 
evident.  Above  3000'’R,  the  chemistry  effect  exceeds  the 
convective  heating  which  eventually  leads  to  model  ignition. 

This  result  is  further  evidence  of  the  validity  of  the 
oxidation  model. 

4.3  Coupling  of  Particle  Erosion  and  Oxidation  Effects 

In  the  preceding  sections,  the  analytical  oxidation 
model  is  shown  to  predict  uncoupled  oxidation  effects  rea- 
sonably accurately.  However,  the  presence  of  particle 
impacts  creates  an  energy  source  which  must  be  considered  in 
evaluating  the  coupling  of  particle  erosion  and  oxidation 
effects. 

A first  order  modeling  approach  was  initiated  due  to 
the  complexity  of  the  particle  impact-oxidation  coupling 
problem  to  attempt  to  understand  the  importance  of  various 
modeling  parameters.  In  conjunction  with  the  one-dimensional 
oxygen  consumption  and  oxide  formation  model  described  in 
Section  4.1,  a 1-D  particle  impact  cratering  and  kinetic 
energy  transfer  model  is  assumed  neglecting  2-D  conduction 
effects.  In  this  regard,  the  model  is  conservative  tending 
to  the  worst  case  situation.  The  steps  in  this  discrete  1-D 
modeling  approach  include; 

1)  Oxygen  dissolution  and  oxide  formation  with 
transient  1-D  conduction  and  diffusion  during 
the  first  particle  obscuration  time. 

! 

2)  Discrete  particle  impact  and  removal  of  a 
1-D  slab  of  oxide/titanium. 

3)  Deposition  of  a fraction  of  the  incident 
particle  kinetic  energy  in  the  titanium  as 

an  instantaneous  distributed  heat  source  over 
a depth  of  the  order  of  the  particle  diameter. 


4)  Transient  1-Ii  conduction  and  oxygen  diffusion 
with  oxide  re-formation  over  the  next  particle 
obscuration  time. 

This  discrete  modeling  approach  results  in  local 
crater  hot  spots  due  to  kinetic  energy  deposition  followed 
by  thermal  diffusion  {and  thus  crater  cooling)  over  the  time 
between  impacts.  When  the  bulk  temperature  reaches  a cri- 
tical level,  a subsequent  impact  causes  a local  ignition  due 
to  the  high  local  crater  temperature  and  thus  chemical  heat 
release  locally.  During  this  latter  phase,  the  solid  phase 
oxidation  and  dissolution  process  becomes  so  large  as  to  be 
gas  phase  diffusion-limited.  Thus,  the  boundary  layer 
transfer  processes  also  become  coupled  to  the  overall  solid 
phase  oxidation  process. 

It  must  be  emphasized  that  the  model  presented  herein 
for  this  highly  coupled  particle  erosion/oxidation  ignition 
process  is  intended  to  identify  the  key  features  involved. 

In  this  respect,  a number  of  adjustable  parameters  have  been 
introduced  relating  to  the  details  of  the  particle  kinetic 
energy  deposition  process  and  to  the  gas  phase  diffusion 
limited  oxidation  rates. 

The  key  parameters  in  this  modeling  include; 

Ppoo  ■ particle  mass  density 

dp  5 particle  diameter 

Vp  H particle  velocity 

6 = titanium  thickness 

Q = cold  wall  convective  heat  transfer  rate 
cw 

H BL  recovery  ethalpy 

< ^ 

Ujj  = depth  of  kinetic  energy  deposition  as 
fraction  of  dp 
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n = fraction  of  particle  kinetic  energy 
transferred  to  thermal  energy 

m^  = augmented  BL  oxygen  mass  flux  normalized 
^ by  unaugmented  value 

Table  4-4  on  the  preceding  page  illustrates  the  represen- 
tative values  chosen  for  this  parametric  study.  The  results 
of  the  parametric  calculations  are  shown  in  Figure  4-17.  No 
local  hot  spot  ignition  is  predicted  except  for  case  3 with 
1mm  particles,  0.1  in.  thick  titanium,  and  an  augmented  BL 
oxygen  mass  flux  by  a factor  of  2.  For  unaugmented  mass 
flux,  the  problem  is  gas  phase  diffusion  controlled  which 
prevents  an  ignition  in  this  case.  Although  case  3 pre- 
dicts, qualitatively,  ignition  for  larger  particles,  cases  4 
and  5 illustrate  that  this  can  be  alleviated  using  thicker 
titanium  material.  These  results  seem  to  indicate  the 
importance  of  the  gas  phase  oxygen  mass  flux  diffusion 
limitations  in  preventing  ignition  at  higher  temperatures. 
The  results  also  indicate  an  increased  concern  for  ignition 
for  higher  dust  densities  and  particle  sizes,  but  ignition 
was  not  predicted  for  conditions  of  interest  below  approxi- 
mately 2000‘’F  bulk  temperatures  which  still  appears  to  be  a 
reasonable  design  criteria. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  conclusions  based  on  this  exploratory  program  to 
evaluate  the  combined  effects  of  simulated  ice/water  and 
dust  environments  on  titanium  surface  chemistry  can  be 
summarized  as  follows: 

(1)  Phase  I test  data  show  v/ater  vapor  environments 
produced  increased  oxidation  rates  by  approximately  100  percent 
for  titanium  compared  to  O2  environments  for  1200-1500°K  but 
similar  rates  for  1700°K. 

Analytical  calculations  suggest  changes  in  H,0  vapor 
oxidation  rates  with  temperature  may  be  explained  based  on 
competing  effects  of  diffusion  coefficient  (increased  H2 
diffusion  coefficient  with  temperature  in  titanium)  and 
solubility  (decreased  H2  solubility  with  increasing  temp- 
erature) . 

(2)  Phase  I test  data  shows  apparent  pressure  effect 
on  oxidation  rate.  Analytical  modeling  where  oxidation  is 
solid-phase  diffusion-controlled  contradicts  this  experi- 
mental result,  x^ssessment  of  literature  data  also  contra- 
dicts this  experimental  result.  The  possibility  of  gas- 
phase  rate-limiting  processes  in  the  Phase  I test  results 
may  be  an  explanation.  However,  the  overall  evidence  from 
the  literature  suggests  no  pressure  effect. 

(3)  Phase  II  combined  H2O  and  dust  data  in  the  J^EDC 
DET  suggest  a dominance  of  the  heat  absorbed  in  vaporizing 
the  H2O  over  any  heat  released  due  to  surface  oxidation. 

The  overall  effect  is  a net  cooling  due  to  the  presence  of 
condensed  H2O.  Application  of  the  analytical  oxidation 
model  to  the  DET  test  data  shows  good  agreement  predicting 
the  on-set  of  significant  oxidation  heat  release  at  approxi- 
mately 2000°F  with  ignition  occurring  at  temperatures  in 
excess  of  2500®F. 
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(4)  Phase  III  tests  comparing  ice  and  glass  single 
impacts  on  titanium  in  oxidizing  environments  showed  ^ 
ignition  for  a bulk  temperature  of  1700°F  ('^^1200°K)  . The 
ice  produced  a smaller  crater  than  the  glass,  but  similar 
oxide  coating  spallation. 

The  overall  conclusion  from  this  exploratory  program 
is  that  nuclear  cloud  ice/water  will  not  significantly 
affect  titanium  oxidation  response.  Design  criteria  de- 
veloped from  extensive  dust  environment  studies  are  still 
valid.  No  ignition  is  expected  if  titanium  bulk  tempera- 
tures remain  below  approximately  2000“F.  Ignition  can  be 
expected  above  approximately  2500°F. 
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